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Abstract We estimate the growth in time of the solutions to a class of non- 
linear fractional differential equations Dq_^_ ( ) = f{t,x) which includes 
Dq_^_ (x — xo) = H{t)x'^ with A e (0, 1) for the case of slowly-decaying coef- 
ficients H. The proof is based on the triple interpolation inequality on the 
real line and the growth estimate reads as x{t) = o{t°'°') when t — )■ +oo for 
l>a>l — a>A>0. Our result can be thought of as a non-integer 
counterpart of the classical Bihari asymptotic integration result for nonlin- 
ear ordinary differential equations. By a carefully designed example we show 
that in some circumstances such an estimate is optimal. 
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1 Introduction 

The fractional calculus deals with the generalization of differentiation and 
integration to non-integer orders. Fractional calculus is as old as the classical 
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differential calculus and it gained more and more importance during the last 
three decades in many fields of science ad engineering [521 EZl ESI E31 EDI ESI 

EIEIIEZIEHIES]. 

The fractional order differential equations play a significant role in model- 
ing the anomalous dynamics of various processes related to complex systems 
in most areas of science and engineering. However, only a few steps have 
been taken toward what may be called a coherent theory of these equations 
in the applied sciences in a manner analogous to the ordinary case. 

We would like to recall that the fractional differential equations have been 
found effective to describing various physical phenomena such as damping 
laws, rheology, diffusion processes [191 HH] ■ Various applications of fractional 
differential equations in the reaction kinetics of proteins, the anomalous elec- 
tron transport in amorphous materials, the dielectrical or mechanical re- 
laxation of polymers, the modeling of glass-forming liquids and others are 
presented in [2S1 HZl ElES] . 

On the theoretical side, several fractional variational principles were for- 
mulated and the corresponding Euler-Lagrange equations consisting of left 
and the right derivatives were analyzed during the last few years [H IH El E2] • 
The fractional derivatives are also the infinitesimal generators of a class of 
translation invariant convolution semigroups which appear universally as at- 
tractors [381 EH] . 

New methods for solving the nonlinear fractional partial differential equa- 
tions were developed recently [2S1 EH] ■ 



Let us consider the initial value problem for a fractional differential equa- 
tion (FIVP) below 



where the nonlinearity f : D = [0, +oo) x M — > M is assumed continuous and 
Xq is some real number. 

The differential operator Dq_^ in problem ([TD is the Riemann-Liouville 
differential operator of order < a < 1, namely 



where r(l — a) = /q °° e~^t~°'dt is the Gamma function. See [2T1 p. 70]. 

The above fractional differential equation represents a wide class of dis- 
sipative nonlinear nonlocal equations, therefore the asymptotic behavior of 
its solutions is an interesting topic of investigation. 



Do\(x-xo)(t) = /(t,x(t)), t>0, 

X(0) = Xq, 



(1) 
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The motivation for inserting the initial datum xq into the differential 
operator comes from the physical origin of such mathematical models and 
the reader can find comprehensive details in this respect in [33l p. 80], [131 
p. 230]. The Riemann-Liouville operator with an inserted datum is called a 
Caputo differential operator — see [9], [2T], p. 91]. 

Assuming that the FIVP has a solution x{t), the formulas r(Q;)r(l — a) = 
. and 

sin Tra 

see [21 p. 196], allow us to rewrite ([T]) via an integration as 

y„ r(l-a)(«-s) 
where t > 0. 

We can regard from now on, by means of ([2]), any solution of FIVP ([T]) 
as a (continuous) solution of the singular integral equation 

1 /"* fis, x(s)) , , , 

and vice versa. This fact is established in a rigorous manner in [13j. See 
also the monograph [211 Theorem 3.24, pp. 199-200] (for 7 = r = in the 
original notation). 

Given the generality of /, we can replace it with r(a)~^/, which means 
that we are interested in analyzing the next integral equation 

xit)=Xo + l y^l'f^Jj s, t>0. (4) 

Our aim here is to investigate the long-time behavior of the solutions to 
a large class of equations (jl]). Noticing that, for a constant / in (jl]), the 
solution reads as 

f 

x(t) = xo + --t" = Oit") as t -> +00, 
a 

we shall compare the growth of solutions to the nonlinear integral equation 
(jH) with the natural test-function t". 

To the best of our knowledge, the literature of fractional differential equa- 
tions is scarce with respect to such enterprises. Most of the attempts regard 
exact solutions to (very) restricted types of fractional differential equations. 



x(s) — Xq 



r(a) 



[S — T 



1 1— a 



dr 



ds = 0, 



(2) 
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E.g., in [T7] the linear fractional differential equation Tq°'Dq_^_ + — 
$(0) =0 has been used to model the complex protein dynamics. The solu- 
tions were expressed using a Mittag-Leffier transcendental function of order 

a, namely <l>{t) = $(0)E„ (^)°), see [m p. 49]. So, ~ t"" as the 
time t increases indefinitely. 

The paper is organized as follows. In section[2]we sketch an asymptotic in- 
tegration theory for nonlinear ordinary differential equations of integer order 
that motivates our investigation. In section |3] we introduce the hypothe- 
ses concerning the functional quantity f{t,x) from ([3]). The main result is 
established in section HI Section |5] contains a detailed presentation of the fun- 
damental class of nonlinear differential equations given by f{t, x) = H{f)x^, 
with a > A > 0. The last section is devoted to an analysis of the sharpness 
of our asymptotic estimate for the continuable solutions to ([3]). 

Let us conclude this section by stating that, even though we have chosen 
to use the simple Riemann-Liouville differential operator inside the nonlinear 
fractional differential equation ([1]) for reasons of clarity in the exposition, the 
technique presented here works equally well for most of the (far more compli- 
cated) differential operators of non-integer order that can be encountered in 
both the mathematics and physics literature. Another point that has to be 
emphasized is that, due to the interpolation-type approach described in this 
paper, we allow the coefficient H to be more flexible (both in the vicinity of 
the initial moment and for the large values of t) than it is usually allowed 
in the recent investigations based on the (weighted) Holder-type inequalities. 

2 The Bihari asymptotic integration theory 

In two seminal papers [3 |S] on the asymptotic integration of nonlinear ordi- 
nary differential equations of integer order, I. Bihari established that, given 
the equation 



and /, a, g are continuous functions, g is monotone non-decreasing and 



x" + f{t,x) = 0, t>to>0, 



where 
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all the continuable, real-valued, solutions x satisfy the conditions 

x(t) 

lim x'{t) = lim = l{x) E M, 

i— >-+oo t^+oo t 

that is, they resemble with a certain degree of accuracy to the lines l{x) ■ t. 

For the case l{x) = 0, of interest for applications in mathematical physics, 
see also the recent contribution [3T] . 

The core of Bihari's approach relies on using his famous integral inequal- 
ity, see the presentation from pQ, to establishing that, if a solution x is 
globally defined in the future, the quantity 

x{t) 

is bounded in [to,+oo). 

This implies further that the improper integral 

/■ + 00 ft 

/ f{s,x{s))ds = lim / f{s,x{s))ds 
J to *^+°° J to 

= x'{to) - lim x'{t) 

t— >+oo 

is convergent, its value being x'(to) — l{x). 

To search for possible answers to the following question — are there 
any non-integer counterparts of the classical Bihari asymptotic integration 
result? — , we have thus (1) to produce a quantity like ^ that will ultimately 
be bounded and further (2) to establish that lim x'{t) is finite. 

t-^+oo 

In the following, the quantity from (1) will be 

x{t) 



B{x,t) 



The part (2), however, cannot be translated into the non-integer case: 
this is a consequence of the non-local nature of the differential operator Dq_^_. 

To answer affirmatively to the preceding question, we have replaced the 
finiteness of lim x'(t), which is equivalent with the finiteness of lim 

with a condition in the same spirit, namely 



lim t^^-'^)" ■ , 



for some a G (0, 1). 
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3 Preliminaries 



Several technical requirements for the nonlinearity / will be presented next. 
Given a G (0, 1), set a G (0, 1) and P2 > 1 such that a > 1 — a and 

(5) 



(1— a)(a + a^) 1 (1 — a)a + 
, 1 - a)a = < — < < a. 



1 + a p2 1 + a 

Introduce further pi,P3 > 1 via the formulas 

1 1 _ 1 + (1 - a)a , 1 + a 

Pi P2 

Notice that 



1 + a 



and 



P3 



aa 



111. 
- + - + - = 1, 

Pi P2 P3 



(6) 



which makes Pi,P2,P3 perfect candidates for the triple interpolation inequal- 
ity on the real line, see e.g. [121 p. 146]. For further use, remark that the 
double inequality regarding p2 from ([5]) is equivalent with 



1 1 - (1 - a)a^ 

1 - a < — < ^ — . 

Pi 1 + a 



(7) 



Consider now the continuous, non-decreasing function g : [0,+oo) — 
[0, +oo) with the following features 

g{0) = 0, giO > when ^ > 0, 

and — take 7^ in — 



lim 



W{u) 



+00, W{u) = / -r—T 
[9[ 







P3 


9 {\xq\ - 







An example of g is that of g{u) = (i+C)« for some C ^ (0) 1) 
a G ^0, Y^j- It is easy to see that 



W{u) 



1-(1+C)a 



PS 



> 



\xo\P-i 



i-(i+C)a 



P3 



, , _ ,,(1+C)a _ U„|{l+C)aP3 

2-[i-(iK)a]P3 . ^ l^ol 

(1 + C)« 

c ■ u^^^'^'^°- when u +00, 



where c = 2-^^-^^^^^"^^-^ ■ 7-^. 

(1+C)a 
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Remark 1 A crucial fact should be emphasized at this point. We shall use 
in sections El El the function g given by g{u) = for a > 1 — a > X. This 
restriction upon A, which can be recast as 1 — A > a, is natural in the above 
computation since we can take (1 + C)a = l — A for a ( small enough. 

We introduce next the class of Bihari-like nonlinearities / by means of 
the inequality 

\f{t,x)\<h{t)g(^-^^^y (9) 

where the function h : [0, +oo) — )■ [0, +oo) is continuous and such that 

t-''""^^'""^^ ■ IW^..(o,) = 0(n whent^+oo. (10) 

Remark 2 The ratio inside g from is designed, following the tradition 
of [SI pp. 277-278], according to our aim of studying the so called long-time 
behavior (that is, for t — )■ +oo) of solutions to ([T]). This is why the choice for 
the argument of g looks sloppy when t is close to 0. We emphasize that it is 
only the quest for clarity and simplicity that has lead us to introducing such 
a peculiar element [t + 1)". A careful inspection of the computations in the 
following will allow the interested reader to adapt the description of / for 
the vicinity of 0. This principle is also valid for the important issue of (any) 
solution a; of ([3]) being or not being allowed to take the values ±oo for t = 0. 
Not to complicate the analysis, we opted for asking x to be continuous as 
well in 0. This limitation yields a;(0) = Xq for all the solutions of ([1]), ([3]). 

Remark 3 Our main interest in this paper is about fractional differential 
equations where the nonlinearity f{t, x) is power-like. This is motivated by 
many applied examples: logistic equations used in population dynamics and 
biochemistry (e.g. polymer growth), Emden- Fowler and Thomas- Fermi type 
of equations (from stellar dynamics where the chaotic expansions of the outer 
space gases might be modeled with fractional differential equations). The 
architecture of the integral W from equation ([S]) has been chosen, however, 
as to give a certain freedom to the nonlinearity. One can have even some 
non-deterministic/non-smooth alterations of the power-like form. 

The restriction (|T0|) can be recast as 

t — 

^fa[i_p,(i_„)] _ 1^ lh{s)f'dsy' < M{t + 1)°, t > 0, (11) 
for some sufficiently large constant M. 
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Remark 4 Notice the inequality of the exponents of t, namely 



Pi 



[1 -pi{l - a)] < a, 



which is equivalent with the second part of the double inequality ([7]). This 
means that the class of coefficients h is quite substantial, including, for in- 
stance, the space L^^(0,+oo). It is a common practice in the asymptotic 
analysis of ordinary differential equations of integer order to study carefully 
the case of slowly- decaying coefficients h, e.g. h{t) ~ t"^ for some e G (0, 1) 
when t — )• +00. This issue seems not to be addressed yet in the literature de- 
voted to asymptotic integration of fractional differential equations. For this 
reason, we hope that our interpolation-type approach will lead to further 
investigations in this area of research. 

4 A Bihari— type result for fractional differ- 
ential equations 

Our main contribution in this work is given in the sequel. 

Theorem 1 Assume that the function f from ^ verifies the conditions 
Then, all the continuable solutions (that is, continuous solutions 
defined in [0, +oo) ) of the FIVP IHj, ^ have the asymptotic behavior 



Remark 5 We have refrained here from discussing the delicate matter of 
the global existence of solutions to This requires the tools of functional 
analysis and its insertion would have lead us too far afield. However, a long 
but straightforward adaptation of the procedures from [1] or [15] will conclude 
that such globally defined solutions of the FIVP ([I]), ([3j) do exist. See also 
Lemma [3] of the present paper and [B] for a significant particular case. 

Proof. We split the demonstration into two steps. 

In step 1, we shall establish the raw asymptotic description of solutions 
to ([1]), ([3]) given by the formula 



x{t) = ait'"'') 



as t ^ -\-oo. 



(12) 



x{t) = o{t") when t — )■ +oo. 



(13) 



To this end, consider x{t) a continuable solution of the integral equation 
(jlj) and set y(t) = x{t) — xq for every t > 0. 
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We have the estimates 



m\ < 



\yis)\ + \xo\ 



ds 



< 



f 


1 


pi 


Jo 







s + 1)° 
LP2{0,t) ■ [z{t)]^ , t>0, 



where 

z{t) = 

Further — recall ffTT 



\y{s)\ + \xo\ 

(s + l)" 



1 P3 



ds. 



< Mi{t + iy ■ z{t), 

with c = [1 - - a)] fi and Mi = c- M. 
We deduce now that 



t > 0, 



(14) 



1 \ P3 



as 



z\t) < [gMt)\ + |xo|)r < 9 (kol + [Mi(t + 1)" ■ 
Fix to > 0. Notice that, when t E [0,to], the preceding inequality reads 

/ 1 \P3 

z\t)<g(\xo\ + [ci-z{t)]p3) , ci = Mi(l + tor, 



and 



/ J_\P3 — '^l- 

g (\xo\ + [ciz{t)]p-s j 
An integration with respect to t leads to 

W(ci2(t)) < cit < cito < Mi{l + toy-^\ t e [0,to]. 
Replacing t with to in ( |T5|) . we get — recall f|T^ — 

W{\y{t)\P-') < W{Mi{t + lTz{t)) < Mi{t + t > 0. 

The inequality (|T6|) can be rewritten as 



P3 



ap3 



+ 



aapz 



< Ml. 



(15) 



(16) 



(17) 
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Notice that the third factor from the second-hand side of ([T7|) is 1, that 



IS 



<Mi, t>0. 



Suppose now, for the sake of contradiction, that y{t) is not o(t") when 
t — +00. This means that there exist an increasing, unbounded from above, 
sequence of positive numbers (t„)„>i and e > such that 



\y{tn)\ 



n> 1. 



(19) 



Since f|T9l) imphes that hm \yitn)\ = +00, we have — recall ([8]) 

71— > + 00 



lim 



lim 



W{u) 



-00. 



According to the inequality (ITHl) . we obtain that lim ^-jj^ = and this 
contradicts (IT^ . 

Step 2. According to the first step, y{t) = o(t'^) for all the large values of 
t. This yields, via the L'Hopital rule, that 



lim 

t^ + OO t 



lim 



Jo 



9 



\y{s)\+\xo\ 

(s+l)« 



P3 



ds 



lim {m±M 



= ^?(o) = o. 

Further, by means of the estimate (fill) , we conclude that 
= o (^t^j = o(t"") whent-^+cx). 



The proof is complete. □ 



5 Applications 

To illustrate this work, let us fix a, A G (0,1) with a > A. There exists 
C € (0, 1) with the property that 

a > (20) 

Set also a = a{Q = j^. 
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The number p2 > 1 satisfying (jSj) can now be introduced. Noticing that 
the function ( i— )■ = 1 — a{() is increasing, the conditions (E]) can be 



i+C 

simphfied using the stronger restriction 

Aa < -■«< — < Aa+ll — A) 

1 + C P2 1 + a 



Aa + a 
1 + a 



(21) 



< 



1 + a 



In fact, let us fix p2 such that 



1 

Xa < — < Xa+ (1 - A) . 

P2 1 + a 



(22) 



Then, we can replace fl20l) . fl2Tl) with the more practical requirements for C 
given by 

A+C A+C 1 
a > -, ■ a < — . 

1 + C 1 + C P2 

Now, the data regarding the numbers a, {pi)i<i<3 can be summarized in 
the following algorithm- like structure: 
Step 1 Take a, A such that 

1 > a > A > 0; 

Step 2 Take p2 such that 



Aa + (1 - A) 



a 



1 + a 



<P2< (Aa) ^ 



Step 3 Take C such that 



Step 4 Set 



Step 5 Set 



< C < mm <^ , \ 

1^ 1 — a ap2 — 1 J 



1 - A 



1 + C 



Pi 



l + (l-a)a 1 
1 + a p2 
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Step 6 Set 



P3 



1 + a 



aa 



Consider the initial value problem for a fractional (general) logistic equa- 
tion with sublinear exponents 



D^+(a; - xo) 
x(0) = Xo, 



H{t)x^ [1 + P{t)x 



where /i G (0, A), the function P : [0, +oo) 
and 



t > 0, 



is continuous and bounded, 



H{t) = {t + iy^h{t), he LP^iO,+oo). 

Then, all the positive-valued, continuable solutions of the initial value 
problem have the asymptotic behavior x{t) = o{t°') when t — )■ +oo, according 
to the first step of Theorem [TJ This lack of sharpness in the asymptotic 
estimate is explained in the next paragraph. 

Notice that 

\H(t)x^[l + P{t)x-^]\ < + ■ (1+ ||P||oo)(a;^ + x'') 

l + x 



< 2{l + \\P\U)h{t) 

< h{t)g 



X 



t > 0, X > 0, 



where = A — /U G (0, 1) and 

g{u) = 2{l + \\P\\^){l + u)\ 

Since this function g does not verify the hypothesis ^'(0) = 0, we can extract 
from the inferences of Theorem [1] only the conclusion of the first step. 

Let us observe also that for a, A and p2 subjected to the preceding re- 
strictions, that isO<A<a<l and p2 confined by (1221) . all the continuable 
solutions of the sublinear fractional differential equation 



Dq+{x - Xo) 
x(0) = Xo, 



H{t)x^, t > 0, 



(23) 



where 



i/(t) = (t+ l)--^/i(t), he LP'iO,+oo), 
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have the asymptotic behavior x{t) = o{t"'°') when t — ?► +00, according to 
Theorem [H 

As detailed in section [3l the asymptotic integration results for these initial 
value problems are manageable illustrations of our technique. Recently, their 
superlinear variant, namely the case when P = 0, a G (|, l) and A > 1 has 
been discussed in the beautiful contribution [15], however, under a quite 
demanding restriction on the coefficient h — for instance, h G i^^(0, +00) 
for all p G (0, 1). Such a restriction excludes all hope for an analysis of the 
long time behavior of solutions in the case of slowly-decaying coefficients. We 
emphasize the fact that, even though most real-life applications of (fractional) 
differential equations involve nonlinear equations, in the research literature 
such enterprises are scarce. In this respect, the contribution [15] is quite 
singular since its methods are mathematically rigorous. On the formal level, 
a solution of type 0{t"^) as t — )■ +00 for a = | and A = 4 has been presented 
in [331 P- 236] for a fractional model of radiative cooling of a semi-infinite 
body. 

6 A direct analysis 

The final section deals with the crucial issue of sharpness for our asymptotic 
estimate in Theorem [TJ We claim that the asymptotic estimate given in this 
paper by means of a Bihari-type approach in conjunction with an interpola- 
tion inequality is almost optimal. To establish the validity of our claim, we 
use a method frequently encountered in the asymptotic integration theory 
of differential equations of integer order, namely we build an "extreme" ex- 
ample and show that, with respect to it, our theorem does the best job ever 
possible. 

Before proceeding with the computations, let us mention that, similarly to 
the case of the classical Bihari investigation [1] , we obtain here the following 
type of error: for an exact solution having the numerical values described by 
0{t'^''^) when t — )■ -|-oo, where e G [eq, 1) for a prescribed £0 > and g > 0, 
the asymptotic behavior obtained using inequalities and averaging reads as 
o{t^) for if: — )■ -|-oo. Since we can set Eq as close to 1 as desired, it is obvious 
that the estimate provided by such a Bihari-type analysis is optimal. 

Consider now the integral equation below 

x{t)=xo + ^^^^^^\_J x{s)]^ds, t>0, (24) 
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where Xq > and a, X E (0, 1) are fixed numbers such that 

1 > Y and a> 2A. (25) 

Here, H : [0, +oo) — )■ [0, +oo) is a continuous function. 

The next lemmas address three fundamental issues regarding the inves- 
tigation of (^^: the sign of x, its uniqueness and its global existence in the 
future. 



Lemma 1 If x : [0,T) — )■ M, where T < +oo, is a continuous solution of 
equation ([^p such that x{0) = xq then 

x(t)>xo, te[0,T). (26) 

Proof. Notice that it is enough to prove that x takes positive values 
everywhere in [0,T). In fact, if we assume that x has a zero in [0,T) then, 
given its continuity, there exists the number to > such that — recall that 
x{0) = xo> — 

x{t)>0 for alH G [0,to), x{to) = 0. 
By taking into account (El]), we get 



x{to) = Xo+ TT— :F(-S)J ds > Xq. 

Jo 1^0 — 

We have reached a contradiction. 

Finally, since x{t) > throughout [0,T), we obtain fl2^ . □ 

Lemma 2 Given the number xo > 0, there exists at most one continuous 
solution X : [0,T) — )■ M, whereT < +oo, of equation [24^ such that x{0) = xq- 

Proof. Assume that there exist two continuous solutions x, y : [0, T) — t- M 
of equation with a;(0) = y{0) = Xq. Then, via Lemma [H we have 

x{t),y{t)>xo for alH e [0,r). 

The mean value theorem, applied to the function A; : +oo) — > (0, +oo) 
with the formula k{^) = yields 



/ 2 \ 

■ {max{x,y} -mm{x,y}y 



< A — -Ix-yl, 
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where ^ lies between x and y. 

The latter estimate implies that 



^ c-^ te[0,T), (27) 

where c = A(2xq "'^)"'^^'^. 

The Lipschitz-type formula (l27|l allows one to apply the iterative proce- 
dure presented e.g. in [33l pp. 127-131] to conclude that the solutions x, y 
coincide. □ 

Lemma 3 Given the number xq > 0, there exists at least one continuous 
solution X : [0, +oo) — )■ M o/ equation [24^ such that x(0) = Xq. 

Proof. The argument is based on a technique developed by the authors 
in [HI Theorem 2]. 

First, set the numbers p,q,L > such that 

11,, . r 1 1 

H — = 1, 1 < p < mm 



p q [a 1 — a 

and 

L > C{a,py for C{a,p) = p^"" • r(l -p(l - a)). 
Second, introduce the functions — here, f{t, x) = H{t)x^ - 

\f{s,xo)\ 



and 



k{t) = l+xo+ [ jl^^^ds, F{t) = c ■ H{t), 
Jo ['^ ~ ^) 

HL{t) = k{t) ■exp(^t + ^ [k{s)F{s)Yds^ , t > 0, 

where c has been defined in the preceding proof. 

Now, consider the complete metric space X = (X^di), where X is the 
set of all elements of C([0, +oo), [xq, +oo)) that behave as 0{HL{t)) when t 
goes to +00 and 

dL[x,y) = s\XY>< r— — \, x,yeX. 

t>o I ^L\t) J 
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Notice that, as in Lemma [H the following inequality holds: 

\f{t,x{t)) - f{t,ym < F{t)\x{t)-y{t)\, t>0,x,yeX. 
Introduce the operator T : X ^ X with the formula 



r{x){t) = XQ + 



f{s,x{s)) 



X eX,t>o. 



Following verbatim the computations from [6J, we conclude that T is a 
contraction of coefficient C{a,p) ■ L^^ which implies, via the Banach Con- 
traction Principle, that our problem has a solution globally defined in the 
future. □ 

Consider now the particular case of equation given by 



H{t) = {t + iy^h{t) = {t + ■ t~~, t > 1, 

where e G (0, 1) is small enough to have (1 — X)a > The existence of e 
will be explained in the next paragraph. As regards the definition of H , this 
function is continued downward to t = in such a way as to keep its sign 
and continuity in [0, +oo). 

It is supposed that the number p2 verifies the double inequality fl22l) . To 
see how to choose e, notice that the inequality below 



1 - \)a > Aa + (1 - A) 



a 



1 + a 



> 



P2 



is equivalent with 



The hypotheses (125|) lead to 



1 > A(2 + a). 



3a 

!>->« 



a 



■ (2 + a) > A • (2 + a) 



a 

Y ~ 2 

so such a number e exists truly. 

Let now x be the solution of equation (l24l) given by Lemmas [H El [3] for 



r = +00 — such a solution exists always 
formula of H, we obtain that 



By taking into account the 



> 



Xo + 




x{s) 




.(s + l)°_ 


Xo + 


f S P3 


Xq 


L it-sy-'^ 


.(s + l)°_ 




ds 1 




I it 


^) t P3 


[t + iy^ 


c . t^'- 


p-i when t — +00 



n A 



ds 



n A 



ds 



{2i 
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where c = — . 

a 

According to our theorem, x{t) = o{t°'°') for t — > +00, where A < 1 — a < 
a. Since we can fix a in this range freely, assume that a is taken such that 



a>2A>l-a>A. 

This choice is in perfect agreement with the hypotheses ([25 
The restriction fl29ll leads to — recall fl22l) — 



(29) 



— > Aa > (1 — A — a)a. 

P2 



Further, we have 



aa > (1- \)a - — > (1 - A)a - ^— 1£ > 0. 

P2 P2 

We would like to evaluate the difference between the exponent aa of t 
in the asymptotic formula of the solution x given by our theorem and the 
exponent (1 — A) a — ^ provided by (l28l) . This reads as 



aa — 



(1- A)a- 



1 +£ 

P2 



l + e 



P2 



+ (a + A - l)a 



<a + e) 



Aa + (1 - A) 



a 



1 + a 



+ (a + A - l)a 



= a ([a + (2 + £)A - 1] + - A)]l 

[ 1 + a ) 

<a{[a+{2 + 6)X - 1] + (1 + e)a} = fia. 



Fix now 770 G (0, 1) and replace the first of hypotheses fl25|l with a stronger 
condition, namely 

7a 

vo>Y- 

Thus, we have 

Q < Tjoa. 

In fact, this inequality is equivalent with 

(1 - 770)0 + (2 + £)A + (1 + £)a < 1. 



Asymptotic integration of fractional DE's 



18 



The latter inequality is valid as a consequence of the estimate below — 
recall the second of hypotheses — 

(1 - r]o)a + (2 + £)A + (1 + e)a < (1 - Vo) + 3A + 2q; 
3a; 

< 1 - r]o + — + 2a < 1. 

The conclusion of this analysis of equation fl2^ is that the solution x, if 
computed numerically, will be at least as big as oo, 
where r] G (0, t/q) and tjq is a prescribed quantity, while the description of x 
given by our theorem reads as o(t"") when t — )■ +oo. In other words, the 
theorem provides a sharp asymptotic estimate for the solutions of ([T]), ([3]) 
and, in some particular cases, this estimate is optimal. 
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